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Abstract 

Ultrarelativistic heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) are thought 
to have produced a state of matter called the quark-gluon plasma, characterized by a very small 
shear-viscosity to entropy-density ratio rj/s, near the lower bound predicted for that quantity 
by AdS/CFT methods. As the produced matter expands and cools, it evolves through a phase 
described by a hadron gas with rapidly increasing rj/s. We calculate rj/s as a function of tem- 
perature in this phase both in and out of chemical equilibrium and find that its value poses a 
challenge for viscous relativistic hydrodynamics, which requires small values of rj/s in order to 
successfully describe the collective flow observables at RHIC. We therefore conclude that the 
origin of the low viscosity matter at the RHIC must be in the partonic phase of the reaction. 



1. Introduction 

Ultrarelativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) are thought 
to have produced a Quark Gluon Plasma (QGP) with the characteristics of a near ideal fluid, with 
large elliptic flow and low shear viscosity to entropy density ratio rj/s [1]. Relativistic viscous 
hydrodynamic calculations require low values of r\js in order to reproduce the RHIC elliptic flow 
{v2> data 12, [2I]. However, current calculations assume a temperature-independent value of rj/s 
throughout the entire evolution of the system. The exact value of rj/s in these calculations is esti- 
mated to be in the range of rj/s « (1 - 3) I An, due to systematic uncertainties, and types of initial 
conditions used [4]. However, it should be strongly noted that the shear viscosity of matter in a 
relativistic heavy ion collision is a time-dependent quantity. While the partonic phase of such a 
collision is expected to have a very low value of rj/s, after hadronization occurs 77/s is expected 
to rapidly increase. Hence, in order to quantitatively constrain the values of rj/s in the deconfined 
phase of a relativistic heavy ion reaction, it is necessary to perform a separate calculation of rj/s 
in the hadronic phase of the reaction, with a systematic analysis of effects of the system evolving 
out of chemical equilibrium in the hadronic phase In these Proceedings, we calculate rj/s 
as a function of temperature in the hadronic phase using the Ultrarelativistic Quantum Molec- 
ular Dynamics (UrQMD) model, both in and out of chemical equilibrium JD]. This suggests a 
potentially strong time dependence of rj/ s in the hadronic phase of the reaction, and hence that 
viscous hydrodynamics practitioners would need to incorporate a temperature dependent rj/s in 
their simulations. 
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2. Model for Hadronic Medium and Viscosity Calculation 



The model for our hadronic medium is UrQMD, a microscopic transport model based upon 
the Boltzmann equation. It contains 55 baryonic species (and their antiparticles) and 32 mesonic 
species, including strange resonances with the most massive baryonic resonance being the iV^sn 
and the most massive mesonic resonance being the </>i9oo (51] • The cross sections used to calculate 
the collision integral in solving the Boltzmann equation are parametrized from fits to experi- 
mental data for processes that are well known and effective symmetries are used to indirectly 
infer cross sections for processes not known well experimentally. The system is initialized with 
particles uniformly distributed in real space, and it evolves through a sequence of binary elastic 
and inelastic processes, and the system is forced into equilibrium by confining the particles in 
our system to a box with periodic boundary conditions. In order to be able to achieve the proper 
chemical equilibrium, we disable multiparticle processes, since it is not possible to implement 
multiparticle collisions using the geometric interpretation of the cross section in UrQMD. Chem- 
ical equilibrium is then verified by comparing hadronic yields to an independent statistical model 
for a hadron resonance gas (SHARE), and kinetic equilibrium is verified by fitting momentum 
distributions to a Boltzmann distribution to extract the temperature. 

In order to extract the shear viscosity of our system, we employ the Kubo formalism. The 
Kubo formalism relates linear transport coefficients to near-equilibrium correlations of dissi- 
pative fluxes and treats dissipative fluxes as perturbations to local thermal equilibrium . The 
Green-Kubo formula for shear viscosity is 



flff(^(0,OK>'(?,f))e q uil, (1) 



where T is the temperature of the system, t is the post-equilibriation time (the above formula 
defines t — as the time the system equilibriates), and n xy is the shear component of the energy 
momentum tensor f7[]. The expression for the energy momentum tensor is 

if = J d 3 p^f(x,p), (2) 

where f(x, p) is the phase space density of the particles in the system. A representative sample of 
the shear viscous correlation function is shown on the left panel of Fig. 1. In order to compuate 
an integral, we assume an exponential ansatz for the correlator, 



Tn 



(jT x H0)7T xy (t)) ocexp\-—\, (3) 



so that the viscosity becomes 



n = ^'(O) 2 )^. (4) 



3. Entropy Considerations 

Caution must be taken when extracting the entropy accurately from a microscopic transport 
model; the system cannot be treated as an ideal gas of massless particles, and specific entropy 
contributions from particles of different masses are not the same. We use the Gibbs formula 
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Figure 1: LEFT: the shear viscous correlator as a function of time for T=67.9 MeV. RIGHT: pressure vs. energy density 
and entropy density versus temperature. 

s = ( e+p zj! BP ' 1 j, and extract the energy density, pressure, baryo-chemical potential, and number 
densities of the relevant chemical species once the system has equilibriated. We verify the cal- 
culation of s via a scaling relation between entropy density and temperature of an equilibriated 

system at fixed volume: s ~ T< 2 < , with c s the speed of sound. The right panel of Fig. 1 shows 
pressure as a function of energy density, and a value of c 2 s = 0.18 extracted from the slope can 
be used to fit the entropy density versus temperature results in that figure. This assures us that 
our calculation is accurate and representative of a hadronic medium of many different hadronic 
species. 

4. Results in Full Chemical and Kinetic Equilibrium 

Using the aforementioned techniques, we present rj/s as a function of temperature in full 
equilibrium for zero baryochemical potential in the left panel of Fig. 2. Also illustrated in that 
figure is the calculation of r//s for chiral pions 101 and 3 flavor perturbative QCD ifioll . The 
minimum value found for rj/s for the equilibrium and zero baryochemical potential case is « 
(0.9-1.0), significantly higher than the KSS bound of rj/s « 0.08. If the minimum value of rj/s 
for hadronic matter in the range of hadronic freeze-out indeed occurs at that value, a viscous 
hydrodynamics fit with a constant value of rj/s to RHIC data may not yield reliable results for 
tj/s in the deconfined phase, since a value of rj/s of at most 0.24 is needed to reproduce RHIC 
elliptic flow data H. 

5. Effect of Chemical Nonequilibrium on i]/s 

However, calculating rj/s in full kinetic and chemical equilibrium, as mentioned earlier, is 
not satisfactory: Although a statistical model analysis of particle yields and ratios at the RHIC 
indicate a chemical freeze-out temperature near r c h e m ~ 160 MeV, hydrodynamic calculations 
indicate that a kinetic freeze-out temperature of T^ n ~ 130 MeV is required to describe the 
momentum distributions of final state hadrons [6]. In order to mimic this situation in our simula- 
tions, we initialize our simulations with off-equilibrium yields and induce a non-unit pion/kaon 
fugacity. The measurement of the viscosity is then made before the system relaxes back into 
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Figure 2: LEFT:^/.v in full chemical and kinetic equilibrium. Also illustrated are the results for 3 flavor pQCD and chiral 
pions. RIGHT: Tj/s at non-unit fugacities and finite baryo-chemical potential. 



chemical equilibrium, and the fugacity is extracted from particle ratios at the time of measure- 
ment. Likewise a finite baryochemical potential can be induced by initializing the system with 
a surplus of baryons relative to antibaryons. Our results for a system with nonunit fugacities 
(red points) and at finite baryochemical potential (green points) are juxtaposed with results for 
tj/s in full chemical and kinetic equilibrium (black points) in the right panel of Figure 2. This 
tells us of a potentially very rich structure for rj/s in a heavy ion collision: as the deconfined 
phase with low rj/s cools and reaches T c , rj/s could possibly rise sharply, then decrease again as 
the hadronic phase evolves out of chemical equilibrium before rj/s rises as it approaches kinetic 
freezeout yjy. Work is now in progress for calculating the bulk viscosity in our model, since that 
is also an effect that is needed for viscous hydrodynamics calculations Ji , IT ]. 
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